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has its own well-defined sky luminance pattern which can be conveniently used to calculate the sky 
radiance and luminance for a given sky patch and the solar irradiance and daylight illuminance on 
inclined surfaces facing various orientations. The crucial issues are whether the skies could be correctly 
identified. This paper presents the work on the classification of the CIE Standard General Skies using 
Keywords: various climatic parameters and indices. Meteorological variables namely luminance distribution for the 
CIE Standard Skies whole sky including zenith luminance, global, direct-beam and sky-diffuse illuminance on a horizontal 
Sky distribution surface and vertical sky illuminance, and horizontal and vertical solar irradiance data are adopted for 
Vertical sky component A F s : 
Climatic parameters analysis. The results demonstrate that there are a number of appropriate climatic parameters for sky 
Sky classification classification and the selection depends on their availability, accuracy and sensitivity. The approaches 
could contribute to the estimation of solar irradiance and daylight illuminance which are essential to the 
renewable and sustainable developments and energy-efficient building designs. 
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1. Introduction 
Solar radiation and natural light data are essential to architects, 
-= Corresponding author. Tel.: +852 34427063: fax: +852 34420427. engineers and scientists for renewable and sustainable develop- 
E-mail address: bcdanny@cityu.edu.hk (D.H.W. Li). ments including active solar energy applications [1,2] passive 
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Nomenclature 


a, b, c,d and e adjustable coefficients 


ay luminous ideal extinction 
By direct-beam daylight illuminance (lux) 
Byr direct-beam illuminance on a vertical surface (lux) 
Dy horizontal sky-diffuse illuminance (lux) 
Dvr sky-diffuse illuminance on a vertical surface (lux) 
E, extraterrestrial illuminance (lux) 
Gy horizontal global illuminance (lux) 
Gyr daylight illuminance incident on a vertical 
surface (lux) 
Ivr sky-diffuse solar irradiance on a vertical plane (W/m?) 
L sky luminance (cd/m?) 
Linea measured sky point luminance (cd/m?) 
lpred predicted sky point luminance in relative form 
Ly zenith luminance (cd/m?) 
My optical air mass 
N the number of reading 
Ryvr ground-reflected illuminance on a vertical 
surface (lux) 
Ts luminous turbidity 


Vis sky illuminance on a vertical surface facing to the 
sun (lux) 
Vos sky illuminance on a vertical surface facing to opposite 


of the sun (lux) 


Z zenith angle of a sky element (rad) 

Zs zenith angle of the sun (rad) 

a altitude angle of a sky element (rad) 

As solar altitude (rad) 

p ground reflectance 

Pwr azimuth angle of the surface normal (rad) 
h azimuth angle of a sky element (rad) 

X scattering angle (degrees) 

Abbreviation 

CIE International Commission on Illumination 
NR normalization ratio 

RMSE root-mean-square error 

SSLD standard sky luminance distribution 

VCR vertical component ratio 

VSC vertical sky component 


energy-efficient building designs [3-6], test reference year databases 
generation [7,8], climatology studies [9-11] and pollution assess- 
ments [12]. An accurate estimation of the available solar radiation 
and daylight is to acquire, not just the total amount of radiation and 
light coming from the sky, but also the radiance and luminance 
patterns over the sky vault [13,14]. The solar radiance and sky 
luminance distributions are influenced by many factors, including 
solar position, turbidity and pollution content of the atmosphere and 
cloud type and amount, which cause numerous interactions of 
sunlight and skylight [15]. Generally, sky conditions can be inter- 
preted into overcast, partly cloudy and clear skies. Sky conditions of 
the same category would have similar solar radiance and sky 
luminance distributions and the corresponding climatic parameters 
would be within certain ranges [16]. Recently, the International 
Commission on Illumination (CIE) has adopted a range of 15 standard 
skies [17] proposed by Kittler et al. [18] as an international standard 
sky model. A number of researchers have reported that the 15 CIE 
standard skies provide a good overall framework for representing the 
actual sky conditions and cover the whole probable spectrum of skies 
found in nature [19-23]. Each CIE General Standard Sky is well- 
defined by sky luminance pattern which is the straightforward 
approach for sky classification [24]. In many locations of the world, 
the basic sky luminance data are not readily available [25]. In 
interpreting sky conditions, meteorological data are initially used as 
weighting factors to show the degree of sky clearness but different 
researchers adopted different climatic indices with various ranges 
[26-30]. The choice depends on the availability, accuracy and 
sensitivity of these meteorological variables [31]. Currently, a number 
of climatic parameters including sky luminance distribution, ratio of 
zenith luminance to horizontal sky-diffuse illuminance (Lz/Dy), ratio 
of horizontal sky-diffuse illuminance to extraterrestrial illuminance 
(Dy/Ey), ratio of horizontal global illuminance to extraterrestrial 
illuminance (Gy/Ey), and vertical sky component (VSC) were pro- 
posed for differentiating the 15 CIE General Standard Skies [32-34]. 
This paper reviews and evaluates the performance of individual sky 
classification processes. Modifications of some sky identification 
procedures are elaborated. A comparative analysis of the perfor- 
mance for various sky classification approaches is reported. Relevant 
climatic elements under individual standard skies are examined and 
implications discussed. 


2. CIE standard skies 


The 15 CIE General Standard Skies are referred to Standard Sky 
Luminance Distribution (SSLD). Each sky standard is the product of 
a gradation function and a scattering indicatrix function [17]. 
The gradation refers to the sky luminance that varies smoothly 
between horizon and zenith. In general, higher sky luminance is at 
the zenith under overcast skies but opposite trend appears 
for clear skies with larger luminance towards the horizon. 
The scattering indicatrix represents the direct sunlight scattering 
in the atmosphere. For a cloudless sky, peak luminance is near to 
the solar position and drops swiftly with the distance from the 
sun [35]. Under overcast conditions, skies are either independent 
of or slight brightening towards the sun. The mathematical 
expressions for the relative sky luminance (ls), relative graduation 
~(Z)/p(0°) and relative scattering indicatrix f(y)/f{Zs) are given as 
follows: 


Et  £We® (1) 
SL fZs)\p(0°) 


(Z) 1+a exp(b/ cos Z) 
0°) 1+a exp b 


(2) 


fo — 1+c[exp(dy)—exp(dz/2)|+e cos? 
f(Zs) 1+clexp(dZs)—exp(dz/2)|+e cos*Zs 


(3) 


where L is the sky luminance in an arbitrary sky element (cd/mz2); 
LZ is the sky luminance at the zenith (cd/m2), Z is the zenith angle 
of a sky element (rad); ZS is the zenith angle of the sun (rad); y is 
the scattering angle (degrees) (i.e. shortest angular distance 
between the sky point and the center of the sun)= 
arccos(cos Zs: cos Z+ sin Zs: sin Z- cos|p—@g|), @ is the azi- 
muth angle of the sky element (rad), AS is the azimuth angle of the 
sun (rad); a, b, c, d and e are adjustable coefficients 

Both gradation and indicatrix functions are of six types. The 
combinations can form 36 skies but only 15 relevant types 
containing five clear, five partly cloudy and five overcast skies 
were adopted to be the standard skies. Table 1 summarizes the 
general features of the 15 CIE Standard Skies. 
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Table 1 
Particulars for the 15 CIE Standard Skies. 


No (code) Gradation Indicatrix 

a b c d e 
1 (11) 4 —0.7 0 -1 0 
2 (12) 4 —0.7 2 —1.5 0.15 
3 (11) 11 -0.8 0 -1 0 
4 (112) 1.1 —0.8 2 —1.5 0.15 
5 (I) 0 -1 0 -1 0 
6 (M2) 0 -1 2 -15 0.15 
7 (M3) 0 -1 5 —2.5 0.3 
8 (4) 0 -1 10 -3 0.45 
9 (IV2) -1 —0.55 2 -15 0.15 
10 (IV3) -1 —0.55 2 -2.5 0.3 
11 (1V4) -1 —0.55 10 -3 0.45 
12 (V4) -1 —0.32 10 -3 0.45 
13 (V5) —1 —0.32 16 -3 0.3 
14 (VI5) -1 —0.15 16 -3 0.3 
15 (VI6) -1 —0.15 24 —2.8 0.15 


3. Climatic variables 


There are a number of daylight parameters for categorizing the 
15 CIE Standard General Skies. The same standard sky type is 
assumed to have the identical well-defined sky luminance pattern 
and the corresponding climatic variables and indices would be 
within certain ranges. Such features can help the identification of 
sky patterns and complicated mathematical expressions to model 
the sky distributions are not required. The characteristics of the 
appropriate climatic parameters for classifying the 15 CIE Standard 
General Skies including L, Lz/Dy, Gy/Ey, Dy/Ey, Ty and VSC are 
discussed as follows: 


3.1. Sky luminance 


Each CIE Standard General Sky is well-defined by sky lumi- 
nance pattern which is the straightforward approach for sky 
identification. Sky luminance distributions particularly for the 
whole sky are seldom available for the sites under consideration 
or exist for very limited periods. The usual instrument 
for measuring sky luminance distribution is by means of a sky 
scanner [35]. There are however, several limitations regarding the 
accuracy of the recorded sky luminance readings. Firstly, the 
sensor head of the scanner rotates in altitude and azimuth to 
measure the luminance at 145 circular sky patches by scanning the 
sky hemisphere. Such a scanning pattern can avoid any double 
counting but it excludes certain regions of sky vault. Secondly, the 
measured data are considered discrete readings rather than 
continuous analytical expressions. Sky luminance between adjoin- 
ing measurement points may vary considerably. Thirdly, each 
scanning period for the 145 points is about 4min and the 
measurements are taken every 10 min. Great changes in sky 
luminance distributions may appear within each measurement. 
Additionally, for “out of range” measurements (points close to the 
sun location under non-overcast skies), an estimation of the sky 
luminance was made from a simple average of the luminance at 
nearby points. Such a conversion however, could introduce data 
distortion. Nonetheless, long-term systematic luminance measure- 
ments for the whole sky can provide sufficient data for analysis. 


3.2. L/Dy 


The ratio Lz/Dy is the original criterion to classify the sky 
luminance patterns into one of the 15 CIE Standard General 
Sky distributions [18]. The integration of the luminance of each 
sky patch over the whole sky vault gives Dy. Zenith is normal to 


Sky description 


CIE standard overcast sky, steep luminance gradation towards zenith, azimuthal uniformity 
Overcast, with steep luminance gradation and slight brightening towards the sun 
Overcast, moderately graded with azimuthal uniformity 

Overcast, moderately graded and slight brightening towards the sun 

Sky of uniform luminance 

Partly cloudy sky, no gradation towards zenith, slight brightening towards the sun 
Partly cloudy sky, no gradation towards zenith, brighter circumsolar region 

Partly cloudy sky, no gradation towards zenith, distinct solar corona 

Partly cloudy, with the obscured sun 

Partly cloudy, with brighter circumsolar region 

White-blue sky with distinct solar corona 

CIE standard clear sky, low luminance turbidity 

CIE standard clear sky, polluted atmosphere 

Cloudless turbid sky with broad solar corona 

White-blue turbid sky with broad solar corona 


the ground surface and Dy would be strongly related to Lz. A good 
agreement (i.e. the zenith luminance is the main component 
contributing to Dy) does exist between these two daylight vari- 
ables [36]. The 15 theoretical Lz/Dy curves can distinguish quite 
well between skies at lower as of 35° or less with clear-cut values 
from 0.1 to 0.41. Sky classification was conducted for high latitude 
regions during winter period (i.e. as less than 35°) [37,38]. 
However, these 15 curves are not parallel and eventually they 
cross each other when as > 35°. This could give ambiguous results 
in sky classification for places of low latitudes [39]. 


3:3. GE, 


The ratio Gy/Ey is often adopted to indicate the relative 
clearness of the atmosphere for sky categorization. In general, 
when the atmosphere is clear, a small fraction of the daylight 
illuminance is scattered, resulting in direct sunlight being pre- 
dominant with a large Gy/Ey reading. Under overcast skies, a large 
portion of daylight illuminance is dispersed, indicating Dy is the 
main component with a small Gy/Ey value. The measurement of Gy 
is simple and straight. The recorded Gy/Ey should be more than 
zero and less than unity. Nonetheless, Gy/Ey may sometimes be 
ineffective to interpret the actual sky conditions. In an extreme 
case, one small cloud in a clear sky may keep the sun obstructed 
by slowly moving across the sky (i.e. small Gy/Ey). On the other 
hand, a small aperture in the clouds in an otherwise cloudy day 
may remain open to the sun for a long period (i.e. large Gy/Ey). 


3.4. Dy/Ey 


As the whole sky is considered, Dy should be an appropriate 
daylight variable to indicate the state of the sky. Basically, the 
measurement of Dy requires a shadow-ring to block the direct 
sunlight. The use of shadow-ring for recording Dy always obstructs 
a significant portion of the sky hemisphere. Correction must be 
applied to the measured Dy to obtain the true readings. Owing to 
the anisotropy of sky-diffuse component with its maximum being 
close to the sun, correction model considering anisotropic effects 
should be employed. High Dy/Ey values often appear in partly 
cloudy skies. For relative low Dy/Ey readings, they may denote 
either overcast or clear sky conditions. Once the three typical sky 
conditions (i.e. overcast, partly cloudy and clear) have been 
distinguished, Dy/Ey can be used to further differentiate the over- 
cast and clear skies [32]. 
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Sky contains air molecules, water vapor, dust and aerosols. The 
attenuation of luminous solar energy through atmosphere indi- 
cates the luminous turbidity (T,,), which is an appropriate index for 
elaborating the daylight climates and sky clearness under cloud- 
less sky conditions. The mathematical expression for T, is given as: 
[18,26]: 


_ In(Ey/By) 
amy 


Ty (4) 


where By is the direct-beam daylight illuminance (lux), ay is the 
luminous ideal extinction=1/9.9+0.043 m, [40], my is the optical 
air mass=1/ sin as +0.50572(as+6.07995)~ 193% [41] 

Under cloudless skies, Ty relates to the atmospheric transmit- 
tance in the direction of the direct sunlight to overall skylight 
efficiency [42]. When there is no Bv, Ty becomes infinite for 
overcast skies. It means that Ty cannot be used to single out 
individual overcast skies. However, it has been shown that Ty is a 
good parameter to differentiate between clean and polluted 
cloudless skies [43]. The typical Ty ranges are between 10 and 15 
for partly cloudy skies and below 6.5 for clear days which can be 
adopted to interpret individual intermediate and cloudless sky 
conditions [32]. 
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Fig. 1. VSC for the five unobstructed overcast skies (CIE Standard General Skies 1-5). 
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Fig. 2. VSC for the five unobstructed partly cloudy skies (CIE Standard General 
Skies 6-10). 


3.6. Vertical diffuse illuminance 


With luminance distribution for the whole sky, daylight illu- 
minance on various tilted surfaces such as vertical planes facing 
different orientations can be estimated by integrating the lumi- 
nance distribution of the sky dome ‘seen’ by the surface [14]. 
Vertical sky component (VSC) is defined as the ratio of the vertical 
diffuse illuminance (Dyr) to the unobstructed horizontal diffuse 
illuminance (DV) available at the same point [44]. The sky-diffuse 
illuminance for the CIE Standard General Skies 1, 3 and 5 can be 
determined using double integrals. However, it is not possible to 
compute the exact sky-diffuse illuminance for other CIE skies 
particularly the non-overcast types. Applying numerical techni- 
ques, Dyr and DV can be expressed as [45,46] 


m n 
Dv= © È Licos a; sin a; 6a; 64; (5) 
j=li=l 


m n 
Dvr = x l X Ly cosa; cos (p— Pyr)iai 59; (6) 
j=li= 
(For O<a<a/2, —a2/2 < $ — nr < T/2;=0 otherwise) 
where æ is the altitude angle of a sky element (rad); and @y, is 
the azimuth angle of the surface normal (rad) 
Horizontal diffuse illuminance is not influenced by the 
azimuthal position of sky luminance patch. The VSC facing a 
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Fig. 3. VSC for the five unobstructed clear skies (CIE Standard General Skies 11-15). 


VCR (Vf8/Vos) 


N 


10 20 30 40 50 60 70 80 90 


° 


Scattering angle (7) 


Fig. 4. Vertical component ratios (VCR) for various y. 
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particular orientation for a given solar position under an individual 
CIE Standard General Sky has its own characteristics and accounts 
for a unique value. The scattering angle (y) is the angular distance 
between a sky-element and the actual location of the solar. The 
sun position is always expressed by two solar angular coordinates 
namely zenith/altitude and azimuth angles. It indicates that a 
given y can stand for many combinations of solar zenith/altitude 
and azimuth angles. For computing y, the location of the sky 
elements is at zero altitude angle and azimuth angle of the 
window normal. Fig. 1 plots the y against the VSC for the Standard 


General Skies 1-5 (the five overcast skies). For sky standards 
1, 3 and 5 (i.e. dark overcast skies with azimuthal uniformity), the 
VSC is independent of y. Constant VSC values for these three sky 
standards range between 40% and 50%. For Standard General Skies 
2 and 4, they are slight brightening toward the sun and the VSC 
changes with y. The VSC varies between 27% and 53% for Standard 
General Sky 2 and from 32 to 64% for Standard General Sky 4. It is 
interesting to note that the two curves do not cross each other. 
Similarly, the VSC at various y for the Standard General Skies 6-10 
(ie. partly cloudy skies) and Standard General Skies 11-15 


Table 2 
The sky number in each of the curve groups. 
Curve group number 1 2 3 4 5 
Sky number 1,3 and5 2, 4,6 and 9 7 and 10 8, 11 and 12 13 and 14 15 
Sky scanner 
For luminance 
distribution L.Lz 
measurements 
a 


illuminance meters 


Photocell(Global/Gv) ©] Gy/y 
; Gv-Dv=Bv—> 
Photocell(Diffuse/Dv) 2 (equation 4) Dv/Ev 
Ns 

Photocell (North) Lz/Dv 
Photocell (East) VSC(based on daylight 

VCR (based illuminance) 

on daylight 
Photocell(south) illuminance) 
Photocell (West) 
Pyranometer(Global) 
Pyranometer(Diffuse) 
Pyranometer(North) N 

VCPbased on VSC(based on 

Pyranometer (East) so.ar irradiance) soalr irradiance) 
Pyranometer(South) 
Pyranometer(West) 


Fig. 5. The measuring mechanisms of the solar and daylight variables used for sky classification. 
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(clear skies) are presented in Figs. 2 and 3, respectively. Under 
non-overcast sky conditions, the VSC depends strongly on the sun 
position. Generally, peak VSC appears at low y of around 10°. Then, 
the VSC drops with increasing y up to about 100°. When y is 
between 100° and 140°, VSC is at its minimum value. Then, VSC 
rises slightly with increasing y. It can be observed that the curves 
intersect with one another at certain solar positions but a 
thorough analysis of Fig. 2 revealed that sky standards 7 and 10 
are quite parallel and do not cross each other. Similar pattern can 
also be found for sky standards 6 and 9. As the cloud amount 
further reduces (i.e. clear skies), the solar corona effect becomes 
more substantial. Larger VSC values are resulted for the General 
Standard Skies 11-15 (i.e. clear skies). At small y less than 45° the 
five curves are separate from each other. However, the curves 
converge and eventually they overlap at y above 50°. Similar 
findings were also reported by Alshaibani [47]. Such features are 
important for sky categorization. 

Under dark overcast skies with azimuthal uniformity (i.e. General 
Standard Skies 1, 3 and 5), the luminance distribution is symmetrical 
about the zenith and changes with the elevation above the horizon. It 
means that the vertical outdoor illuminance would be the same for 
all orientations. For non-overcast skies, very large sky luminance 
would occur near to the sun dropping rapidly with the angular 
distance from the solar, particularly for cloudless skies. As pointed out 
by Alshaibani [48], the ratio of the sky illuminance on a vertical 
surface facing the sun to that facing the opposite orientation (Vg/Vos) 
for instance the vertical component facing south to that facing north 
can show some clear distinctions between groups of certain skies. 
Fig. 4 presents such vertical component ratios (VCR) under various y 
which only refer to the vertical surface facing the sun (surface with 
the smallest y) between 0° and 90°. Totally, six distinguish indicatrix 
categories with the V,/V,; ratios ranging from 1 to over 8 can be 
formed. Each category may contain 1-4 sky types that are summar- 
ized in Table 2 [48]. It indicates that ratios of vertical sky illuminance 
and VSC can be applied for the sky categorization. Fig. 5 presents the 
measuring mechanisms of the solar and daylight variables used for 
sky classification. 


4. Sky classification approaches 


Previously, several approaches were used for classifying the 15 
CIE Standard General Skies [20,32,34]. The sky classification 
procedures with some modifications are elaborated as follows: 


4.1. Sky luminance (Approach I) 


Sky luminance distributions for the whole sky vault perhaps 
are the most appropriate data for sky classification. To identify the 
set of standard skies, the luminance distributions of individual sky 
standards were modeled using Eqs. 1-3 and compared with the 
scanned sky luminance readings. The modeled sky luminance 
should be initially normalized to the horizontal diffuse illumi- 
nance by multiplying all the luminance values with the normal- 
ization ratio (NR) as 


Ylinea COS @ sin a da dh 


NR= z 
Ylprea COS @ sin a@ da dh 


(7) 


Table 3 


where Lmea is the measured sky point luminance (cd/m7); Iprea is the 
predicted sky point luminance in relative form (dimensionless). 

Once the integrated diffuse illuminance and Ipreq at each sky 
point has been obtained, the predicted sky luminance (Lpreq) which 
is the product of NR and Ipreq can be easily computed. An 
alternative would be to divide all sky luminance readings by the 
corresponding Lz, but this can cause huge measuring error when 
the sun is near to the zenith [24]. For low-latitude climates (e.g. 
Hong Kong) when the sun is frequently within a small angular 
distance from zenith, a normalization with respect to the Dy 
would, therefore, be more suitable. Once normalized, the perfor- 
mance of each standard sky luminance model was assessed by the 
root-mean-square error (RMSE): 


1 L d Lea 2 
RMSE = pie ) 8 
a(t (8) 


where N is the number of readings (dimensionless). 

The standard sky interpreted is the one with the lowest RMSE. 
The procedures are repeated for every measured luminance scan. 
Ultimately, the frequency of the identified 15 CIE General Standard 
Skies and the corresponding RMSE can be determined. 


4.2. Appropriate daylight variables (Approach II) 


The original criterion for classifying the sky luminance patterns 
into one of the 15 Standard General Skies is the Lz/Dy which can 
characterize the momentary sky brightness. However, the Lz/Dy 
theoretical curves for the 15 Standard skies are not parallel and 
they intersect with each other at the aS of higher than 35°. The 
global daylight components (Gy) of quite straight measurement 
procedures are appropriate data for classifying daylight climates. 
Direct-beam horizontal daylight (By) can be obtained by subtract- 
ing the corrected Dy from the corresponding Gy. The 15 CIE 
General Standard Skies can be categorized by such four horizontal 
daylight parameters. As proposed by Kittler and Darula [26], the 
pair Lz/Dy—Gy/Ey is an effective hybrid daylight-variable to iden- 
tify the three typical sky conditions (i.e. overcast, partly cloudy and 
clear). There are no clear-cut values for Lz/Dy—Gy/Ey to represent 
different sky types. The Lz/Dy—Gy/Ey ranges for the categorization 
of the three typical skies using Hong Kong measured data were 
computed [43] and are shown in Table 3. To further distinguish 
between the five sky types under each of the typical skies 
(ie. overcast, partly cloudy and clear), we proposed the criteria 
for the various parameters including L7/Dy, Gy/Ey, Dy/Ey and Ty and 
are summarized in Fig. 6 [32]. Two sky classification procedures 
using these daylight variables are adopted. When the data are 
measured at aS of 35° or below, the theoretical L7/Dy curves are 
used for categorizing all the 15 CIE General Standard Skies. The 
criteria for classifying standard skies are universal. For the data 
recorded at aS more than 35°, the procedures presented in Table 3 
and Fig. 6 are applied [49]. The necessary data are the basic 
daylight parameters. Lz is the 145th measurement point for the sky 
scanning. The global and diffuse outdoor illuminance on a hor- 
izontal surface can be used to compute the required daylight 
variables. 


The criteria for the classification of the 3 typical sky conditions based on the hybrid daylight parameter Lz/Dy—Gy/Ey. 


Sky conditions 


Overcast 
Partly cloudy 
Clear 


Criteria 


Lz/Dy > 0.3 and Gy/Ey < 0.3 
0.17 < Lz/|Dy < 0.3 and Gy|Ey< 0.3 or Lz/Dy > 0.17 and 0.3 < Gy/Ey < 0.5 
Lz|Dy < 0.17 or Lz/Dy > 0.17 and Gy/Ey > 0.5 
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Overcast sky dataset 


Sky standards 1,2, 3 & 4 


Lz/Dv > 0.32 


Sky standard 5 
Lz/Dv < 0.32 


Sky standards 1 & 3 


Sky standards 2 & 4 


Dv/Ev < 0.18 Dv/Ev > 0.18 
Sky standard 1 Sky standard 3 Sky standard 2 Sky standard 4 
Lz/Dv > 0.38 Lz/Dv < 0.38 Theoretical Lz/Dv Theoretical Lz/Dv 
curve [18] curve [18] 


Partly cloudy sky dataset 


Sky standards 7, 8, 9 & 10 


With direct 


component 


Sky standard 6 
Without direct component 


Sky standards 7 & 9 


Sky standards 8 & 10 


Tv2 12 


Sky standard 7 
Theoretical Lz/Dv 
curve [18] 


Sky standard 9 
Theoretical Lz/Dv 
curve [18] 


Tv<12 


Sky standard 8 
Theoretical Lz/Dv 
curve [18] 


Sky standard 10 
Theoretical Lz/Dv 
curve [18] 


Clear sky dataset 


Tv 


Sky standards 11, 13, 14 & 15 


>4 


—= 


Sky standard 12 
Tv <4 


Sky standards 11, 13 & 14 
Dv/Ev < 0.3 


Sky standard 15 
Dv/Ev > 0.3 


Sky standards 11& 13 
Theoretical Lz/Dv 


Sky standard 14 


curve 


Theoretical Lz/Dv 
18] curve [18] 


Sky standard 11 
Tv <5.5 


Sky standard 13 
Tv>5.5 


Fig. 6. The criteria for the classification of the 15 CIE Standard General Skies using the appropriate daylight indices. 
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4.3. Vertical sky illuminance (Approach III)/irradiance(Approach IV) 


For sky classification, the VCR (Vg/Vos ratios) from the measured 
vertical illuminance facing the four cardinal orientations (i.e. N, E, 
S and W) were determined and categorized into one of the six 
indicatrix groups (i.e. Fig. 4). Each indicatrix group contains only 
one to four sky types. Apart from indicatrix group 5, the VSCs for 
the sky types within a particular indicatrix group are parallel for 
all y. It means that the sky type can further be identified by using 
VSC. For the indicatrix group 5 (Skies 13 and 14), only the VSC with 
the smallest y (usually less than 45°) was used for sky classifica- 
tion. Daylight is the visible portion of solar radiation and both 
solar radiation and daylight illuminance have similar character- 
istics in nature. The same sky classification procedures can be 
applied to solar irradiance data which are more available than 
daylight illuminance [48] (Approach IV). The methodology is 
universal and only the measurements of the horizontal and two 
vertical sky-diffuse components are sufficient to identify the 
standard skies. 


5. Sky classification and performance evaluation 


The measurements of solar irradiance, daylight illuminance and 
sky luminance were made by a measuring station located at the 
City University of Hong Kong. All sensors were installed on the 
rooftop in a position relatively free from any external obstructions, 
and readily accessible for inspection, general cleaning and main- 
tenance. The data collection started just before sunrise and 
finished after sunset every day. All data were measured approxi- 
mately simultaneously in true solar time which facilitates the 
computations regarding solar geometry and the subsequent com- 
parison of the data at other locations. Six pyranometers were used 
to carry out the measurements of global and diffuse irradiance on 
a horizontal surface and the global irradiance on the vertical 
planes facing the four cardinal orientations (i.e. north, east, south 
and west). The irradiance data were captured simultaneously 
twice per second, averaged over 10-min intervals and transmitted 
to the computers for storage. The daylight illuminance data 
namely the global and diffuse components on a horizontal surface, 
and vertical global illuminance facing north, east, south and west 
were measured by six illuminance meters (T-10M) manufactured 
and calibrated by Minolta of Japan. Similar record procedures as 
solar irradiance were adopted for daylight illuminance measure- 
ments. The sky luminance distributions are scanned using a sky 
scanner (EKO MS 300LR). It records the luminance at 145 points of 
sky by scanning the sky dome. The full view angle of the scanner is 
11° allowing each sky patch to be treated as a point source with 
negligible error [50]. The important parts of the sky scanner are 
housed in a weatherproof casing for continuous outdoor opera- 
tion. The data from the scanner are recorded on a microcomputer 
placed inside the laboratory on the top floor. To safeguard the 
sensor, the scanner does not record luminance data larger than 
35 kcd/m? by using an automatic shutter. Details of the solar 
irradiance, daylight illuminance and sky luminance/radiance mea- 
surements can be found in our previous studies [20,35]. To 
eliminate spurious data and erroneous measurements due to the 
cosine effect of the sensors, proper quality-control tests were 
adopted. Details of the test procedures can be referred to our past 
work [31]. The 10-min data measured in 2004 (January to 
December) were employed for the analysis and the quality control 
tests were conducted. It is inevitable that there were some periods 
of missing data for various reasons such as instrumentation 
malfunction, power failure and senor calibration. Considerable 
efforts were made to obtain a continuous record of data. After 
the quality-control tests, around 15,000 data were retained for 


subsequent analysis. The six pyranometers, six illuminance meters 
and the sky scanner form a systematic measurement to collect the 
climatic variables for categorizing the standard skies. Through the 
stringent quality-controls tests, it indicates that the data recorded 
in Hong Kong have the representativeness for skies classification. 


5.1. Sky classification 


Accordingly, the categorizations of the 15 standard skies using 
the four methods were determined. Fig. 7 presents the results 
using the best-fitting sky luminance approach (Approach I). Large 
variations can be observed for individual sky types and a thorough 
analysis of the figure revealed that the frequency of overcast (Sky 
nos. 1-5) and clear (Sky nos. 11-15) skies was quite close. Totally, 


Frequency of occurrence, % 


123 4 5 6 7 8 9 10 11 12 13 1415 
Sky standard 


Fig. 7. Frequency of occurrence of the 15 CIE skies classified by sky luminance data 
(Approach I). 
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Fig. 8. Frequency of occurrence of the 15 CIE skies classified by appropriate 
daylight variables (Approach II). 


Frequency of occurrence, % 


1 2 3 4 5 6 7 8 9 
Sky standard 


10 11 12 13 14 15 


Fig. 9. Frequency of occurrence of the 15 CIE skies classified by VCR and VSC using 
daylight illuminance data (Approach III). 
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about 75% of the skies were categorized into overcast and clear 
conditions. The partly cloudy skies (Sky nos. 6-10) represent the 
remaining 25%. Sky nos. 1 and 3 dominate the overcast sky. Sky 
nos. 7 and 8 are the main sky patterns for the partly cloudy 
condition. For clear sky condition, it is matched by Sky no. 13. The 
findings are in good agreements with our previous work using 
data recorded between 1999 and 2005 [32]. Fig. 8 displays the 
frequency of occurrence of the 15 standard skies using the 
appropriate daylight variables (Approach II). The distribution is 
quite different. Based on the classification method, more than 50% 
skies were identified as clear state. Sky nos. 1, 7, 13 and 14 are the 
four most frequent standards representing 12.4, 10.6, 18.6 and 
13.3%, respectively. Likewise, the frequency of occurrence of the 15 
standard skies according to the VCR and VSC based on daylight 
illuminance (Approach III) and solar irradiance (Approach IV) are 
respectively plotted in Figs. 9 and 10. It can be seen that the two 
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Fig. 10. Frequency of occurrence of the 15 CIE skies classified by VCR and VSC using 
solar irradiance data (Approach IV). 


Table 4 


figures depict similar patterns. This is not surprising, given that 
daylight is the visible portion of solar radiation. Generally, the two 
figures show that the sky standards are quite evenly distributed in 
terms of the three typical sky types (i.e. overcast, partly cloudy and 
clear) under such categorizations. The majority of the overcast 
conditions are Sky nos. 1 and 2. Referring to partly cloudy 
conditions, Sky nos. 7 and 8 are the main sky type. For clear sky 
condition, Sky no. 13 has the largest frequency. Among all the sky 
classification approaches adopted, Sky no. 13 can be considered 
the most representative clear sky type from the CIE Standard Skies 
in Hong Kong ranging between 15.4 and 18.6%. This supports the 
argument that Hong Kong is currently infamous for its poor air 
quality and polluted conditions. All findings show that in Hong 
Kong Sky nos. 1 and 13 are respectively the main sky patterns for 
overcast and clear sky conditions which are consistent with 
previous studies [8,20,27,32,34,44,49,51 ]. 


5.2. Performance evaluation 


For the development of active solar applications, passive solar 
building designs and daylighting schemes, vertical global values 
are always required. To evaluate whether the skies are correctly 
identified, the global components for the four (i.e. N, E, S and W) 
cardinal vertical surfaces were modeled based on the sky types 
classified by the various approaches. The daylight illuminance 
incident on a vertical surface (Gyr) is the sum of direct-beam (Byr), 
sky-diffuse (Dyr) and ground-reflected (Ryr) components. Mathe- 
matically, it can be expressed as 


Gyr = Byr +Dyr +Rvr (9) 


Calculation of Byr is quite simple, given the position of the sun and 
orientation of the surface. For estimating the Ryz, it is assumed that 


Summary of MBE and %RMSE for vertical daylight illuminance based on the skies classified by Approach I. 


Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Year 
North 
% MBE —12 —7.8 —83 —10.2 9.8 52 10.2 Ia 18.5 17 18.8 21.1 12.4 
% RMSE 13.3 13 14.9 17.1 16.7 13.2 17.4 22.8 219 19.3 20.3 22.1 19.7 
East 
% MBE 3.8 32 0.4 4.2 3:5 4.9 4.6 1.8 —0.7 05 —5.2 -8 1 
% RMSE 11.8 11.7 12.4 15 15.2 14.6 16.9 14.3 12.7 12 14.4 13.3 14 
South 
% MBE 1.4 13 3.2 75 14.9 18.1 19.7 14.8 -5.3 —0.8 1.6 —0.8 -6.3 
% RMSE 5.4 9.3 11 153 19 20.3 22.9 20 13.2 8.3 8.6 8.3 15.8 
West 
% MBE 18.7 10.9 73 2.6 0.4 2 24 2.6 6 10.7 16.8 22.1 —6 
% RMSE 23.8 18.5 16.2 12.2 12.3 12.4 12.6 13.2 15.5 16.6 21.8 25.3 15.6 
Table 5 
Summary of %MBE and %RMSE for vertical daylight illuminance based on the skies classified by Approach II. 
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Year 
North 
% MBE 7.7 6.8 7.6 8.9 6.7 1.8 8.2 13.7 14.5 12.2 16 —19 —9.7 
% RMSE 13.3 15.8 16.7 175 16 12.9 17.4 20.5 20.7 171 19.3 21.6 19.2 
East 
% MBE —0.5 —0.2 3.8 8.6 7 8.2 7.4 55 13 2.2 —3.5 —8.2 3.9 
% RMSE 12.4 14.7 19.3 23.7 21:2 19.9 23.7 20.2 14.2 12.8 15 13.7 18.8 
South 
% MBE 2.9 3.5 3.2 3 11.8 14.9 18.4 10.9 -13 13 3.8 —0.5 —31 
% RMSE 5.7 13.8 15.1 16.5 18.3 18.3 23 18.4 15 9.5 11.3 11 16.3 
West 
% MBE 18.1 8.2 43 11 47 5.9 0.6 1.6 1.3 8.4 15.4 —21.8 -2.8 
% RMSE 24 17.6 16.2 17.7 18.2 18.1 18.4 16.9 16.6 15.4 20.9 25.3 17.9 
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vertical surface receives half of the total illuminance reflected 
isotropically from the ground and ignoring the inter-reflection 
between the vertical and the ground surfaces which is a second- 
order effect of lesser importance. Hence, Byr and Ryr can be 
expressed as follows: 


Byr = (By / sin as)cos as cos (@s-— Pn) (10) 


Ryr = pGy/2 1) 


where p is the ground reflectance (dimensionless) 

For determining Ryr, a common average p value of 0.2 is 
applied. Likewise, the same calculation procedures can be 
employed to determine the sky-diffuse solar irradiance on a 
vertical plane (Ivr). The performance was assessed by two widely 
used statistics, mean-bias-error (MBE), a measure of the under- 
lying trend, and root-mean-square-error (RMSE), an indication of 
data scattering. 

Table 4 presents the MBE and RMSE results for the vertical 
daylight illuminance facing the four principle orientations in all 12 
months and the whole year using the Approach | (i.e. sky 
luminance distributions). The daylight illuminance data for 
north-facing surface are under-estimated for all months with an 
annual MBE of — 12.4%. For other vertical planes, the MBEs vary 
between an under-estimation of 22.1% in December for west- 
facing surface to an over-estimation of 4.9% in July for east-facing 
surface. The RMSEs range from 5.4% in January to 22.9% in July 
appearing both at south-facing surface. The annual RMSEs for the 
four vertical surfaces are below 20%. Likewise, the MBEs and 
RMSEs for the vertical daylight illuminance based on Approach II 
(i.e. appropriate daylight variables) are summarized in Table 5. 
Again, daylight illuminance data are under-estimated in all 
months for north-facing surface. For the whole year, the MBEs 


Table 6 


are —9.7%, 3.9%, —3.1% and —2.8% for north, east-, south- and 
west-facing planes. The peak monthly RMSE of 25.3% appears in 
December for west-facing surface. The annual RMSEs do not differ 
in great deal, ranging from 16.3 and 19.2%. The MBEs and RMSEs 
for the vertical daylight illuminance using the Approach Ill 
(sky classified by VCR and VSC using daylight illuminance data) 
and Approach IV (sky classification by VCR and VSC using solar 
irradiance data) are listed respectively in Tables 6 and 7. In general, 
they have similar annual RMSE pattern. The RMSEs for the 
whole year based on Approach III are 14.6, 20.6, 14.1 and 20.9% 
for north-, east-, south- and west-facing surfaces, respectively. 
Using Approach IV, the annual RMSEs are 0.2-4.2% more than 
those based on Approach III. 

In view of the fact that solar irradiance and daylight illuminance 
have similar characteristics in nature, the methods for daylight 
illuminance prediction can also be employed in solar irradiance 
calculations [52]. Similarly, the monthly and annually MBEs and 
RMSEs for the vertical solar irradiance facing the four cardinal 
directions based on the four approaches were computed and are 
summarized in Tables 8-11. Generally, the four approaches have quite 
good performance in terms of the annual RMSEs. The maximum 
annual RMSE is not more than 21%. The results also show that the 
Approach III (Table 10) predicts the solar irradiance the best with the 
lowest annual RMSE of 9.7% appearing at east-facing surface. 


6. Conclusions 


Various approaches for sky classification were elaborated and 
reviewed. The strengths and limitations of the meteorological para- 
meters including luminance distribution for the whole sky including 
zenith luminance, global, direct-beam and sky-diffuse illuminance on 


Summary of MBE and %RMSE for vertical daylight illuminance based on the skies classified by Approach III. 


Jan Feb Mar Apr May Jun 
North 
% MBE 11 2 11 31 2.9 0.1 
% RMSE 6.6 9 9.4 12.8 16.5 15 
East 
% MBE 9 2 6.3 9.3 6.9 8.5 
% RMSE 16.4 15 152 17.8 14.9 18.2 
South 
% MBE 9.9 4.4 3.6 2 9 6 
% RMSE 11.8 14.9 9.8 12.4 14.3 8.5 
West 
% MBE 35 —2.5 4.1 8.9 8.9 8.8 
% RMSE 26.8 18.5 17.5 20.4 17.7 18.8 
Table 7 


Jul Aug Sep Oct Nov Dec Year 
3.3 8.7 6 5 6.2 3.9 45 
18.9 20 11.7 9.4 9.6 6.3 14.6 
6.2 7 7A 10.8 42 6.2 7 
17.9 24 23.5 28.1 24.4 22.7 20.6 
7A 8.3 0.4 6.1 9.6 9.8 —0.9 
10.9 17.3 15.5 14.8 20.5 13.7 141 
52 6 46 —0.9 —10.2 —9.9 2.7 
20.4 24 23.4 23.3 23.4 17.2 20.9 


Summary of MBE and %RMSE for vertical daylight illuminance based on the skies classified by Approach IV. 


Jan Feb Mar Apr May Jun 

North 

% MBE 25.8 6.5 31 0.5 —0.2 2.4 

% RMSE 31.3 17.5 15.6 11.9 11.1 15.7 
East 

% MBE 21.1 TA 9.2 11.9 6.6 11.6 

% RMSE 25 22.5 26.6 26 19.8 28.4 
South 

% MBE 31 3 6.2 23 3.7 6.1 

% RMSE 10.9 19.8 12.2 11.8 9.8 13.7 
West 

% MBE —2.7 —3.4 1.5 5.9 7.1 6.5 

% RMSE 10.3 27 22.6 26.7 26.9 23.2 


Jul Aug Sep Oct Nov Dec Year 
0.7 48 2.7 —0.6 —2.6 -15 0.1 
14.9 16.2 15.5 14.8 13.2 10.5 14.8 
10.4 47 0.4 31 —3.2 —4.3 6 
25.6 22.6 21.6 19.6 19.2 17.2 23.4 
6.8 3.9 0 —3.8 —8.1 —9.5 —2.5 
13.6 14.1 16.2 13.4 28.3 27 15.9 
2.6 3.4 11 -111 —17.6 —19.1 -1.2 


23.2 26 23.8 23.1 29.8 26.5 25.1 
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Table 8 
Summary of MBE and “RMSE for vertical solar irradiance based on the skies classified by Approach I. 


Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov 
North 
% MBE —11.2 il Be 11 2.2 5.7 56 0.7 5.8 aa 13 
% RMSE 14.6 13.6 18.5 15 15.7 17.4 18 16.1 16.4 14.7 16.4 
East 
% MBE -9 3.2 3.6 5.4 3 4.6 4.7 2.9 1.9 -3 —8.6 
% RMSE 13.7 12.3 12.8 17.6 17.7 19.6 21.7 21.7 20 11.2 15.1 
South 
% MBE —11 41 1.1 8 2.9 35 4.6 55 2.7 4.7 —3.1 
% RMSE 9.4 22.6 18.3 17.8 16.6 16.6 171 19.7 20.6 12.2 13.7 
West 
% MBE —143 —0.2 —0.9 2.5 29 41 2.1 15 43 8.1 13.3 
% RMSE 19.7 15.2 15.2 19.8 17.5 17.1 18.6 18.6 15.2 11.9 16.5 
Table 9 
Summary of %MBE and %RMSE for vertical solar irradiance based on the skies classified by Approach II. 
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov 
North 
% MBE —6.1 2.2 43 1:7 4.7 8.2 7A 22 14 2.7 10.1 
% RMSE 14.1 15.8 22.1 17.4 17.8 18.7 19.6 17 18 14.3 16 
East 
% MBE —5.5 53 5.2 8.5 5.3 6.6 5.9 4.7 4.4 0.3 —6.4 
% RMSE 12.9 16.6 16.8 21.9 20.7 20.6 22.2 22 19.4 12.5 14.8 
South 
% MBE —0.1 7.8 1.2 51 0.2 0.5 3.2 3 0.3 —2.6 —0.7 
% RMSE 8.6 22.1 15.9 18.8 19.5 18.5 19.5 20 20.3 TL7 13.1 
West 
% MBE —124 1.6 03 4 T3 6.7 3.9 44 0.6 —4.8 —11.1 
% RMSE 18 16 16.4 21.6 23.4 19.2 19.5 20.8 17 115 15.8 
Table 10 
Summary of MBE and %RMSE for vertical solar irradiance based on the skies classified by Approach III. 
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov 
North 
% MBE 1 2 1 3.4 4.1 32 55 9.4 8.1 Ta -8.3 
% RMSE 6.3 8.6 9 12.6 12.6 13.9 18 19.9 15:6 13.5 12.9 
East 
% MBE 2.5 —2.2 0.1 1.6 0.4 1.4 0.2 1.6 2.6 2 -6.3 
% RMSE 6.3 11.2 6.2 6.7 9.1 92 9.2 11.8 10.2 9.9 12.6 
South 
% MBE 11 21 2.6 43 7.4 9.9 10.5 9.8 6.5 5.3 —7.2 
% RMSE 5.1 10.9 8.2 11.5 11.1 13.8 15.9 16.7 14.6 12.2 13.4 
West 
% MBE 8.3 6.5 2.9 0.5 0.4 0.8 3 3.6 5.3 9.6 16.1 
% RMSE 10.5 14 7.6 8 6.8 8.3 11.7 13:9 13.3 15.9 22.4 
Table 11 
Summary of MBE and %RMSE for vertical solar irradiance based on the skies classified by Approach IV. 
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov 
North 
% MBE 2.7 75 8.6 6.4 153 9.4 10.2 4.8 4.2 3 —4.1 
% RMSE 9.2 13.1 15.1 15 16 18 18.3 16 16.7 13 12.6 
East 
% MBE 541 9.7 8.1 9.4 16.5 8.8 7.2 9.9 137 8.5 25 
% RMSE 10.6 20.7 13.7 17.4 22.3 20.9 22.6 23.8 28.9 17.6 14.7 
South 
% MBE 13.1 11.3 6.4 —0.3 10.4 3.2 3.2 27 43 5.2 5.5 
% RMSE 17.1 21 13.3 13.1 14 13.9 14.8 1.5 16.8 13.2 14 
West 
% MBE -3 3.8 5:5 10.6 18.4 10.3 10.7 9 4.7 11 —6.1 
% RMSE 14.3 16.2 13.9 22.5 21.5 20.2 20.8 18.6 16.9 14.6 13.7 
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a horizontal surface and vertical sky illuminance, and horizontal and 
vertical solar irradiance data for sky categorization were elaborated. 
Totally, four sky classification approaches with some modifications 
were proposed. Ten-minute solar irradiance and daylight illuminance 
data measured in 2004 were used for the study. The frequency of 
occurrence of the complete set of the 15 Standard General Skies was 
determined accordingly. Sky nos. 1 and 13 dominated respectively the 
overcast and clear conditions. For partly cloudy conditions, sky no. 8 is 
the main sky pattern based on Approaches I, III and IV. To evaluate 
whether the skies are correctly identified, the global solar irradiance 
and daylight illuminance for the four cardinal vertical surfaces were 
modeled based on the classified sky types. The performance was 
assessed in terms of mean-bias-error (MBE) and root-mean-square- 
error (RMSE). In general, the predicated vertical global components 
were in good agreements with those measured values. The peak 
annual RMSEs were 20.7 and 25.1% for solar irradiance and daylight 
illuminance, respectively. The findings support that once the skies 
have been correctly classified, good quality solar irradiance and 
daylight illuminance data on titled surfaces facing different directions 
can be estimated for subsequent analysis. 
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